Driven collective instabilities in magneto-optical traps: a fluid-dynamical approach 
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We present a theoretical model to describe an instability mechanism in ultra-cold gases, where 
long-range interactions are taken into account. Focusing on the nonlinear coupling between the col- 
lective (plasma-like) and the center-of-mass modes, we show that the resulting dynamics is governed 
by a parametric equation of the generalized Mathieu type and compute the corresponding stability 
chart. We apply our model to typical ranges of magneto-optical traps (MOT) parameters and find 
a good agreement with previous experimental observations. 



The celebrated method of magneto-optical cooling of 
atoms [T] allowed in recent years the exploitation of new 
and very interesting phenomena in atomic physics, and 
in particular the creation of Bose-Einstein condensates 
[2]- Cold and ultracold gases produced in a magneto- 
optical trap (MOT) provide an unique medium to explore 
the connections between different areas of physics. Due 
to absorption and re-scattering of light j^^, collective 
phenomena of cold gases confined in a MOT became re- 
cently an issue of major interest. Under such conditions, 
the atomic gas in a MOT can be regarded as a fluid with 
a tunable effective charge, which opens a new area of 
plasma physics [BHH]- A very interesting manifestation 
of the latter is the work on the Coulomb explosions of 
molasses, performed by Pruvost et al [7]. Moreover, in 
a recent work, we have proposed the existence of hybrid 
waves and Tonks-Dattner oscillations in the confined gas 
of ultra-cold atoms, very similar but not identical to those 
of a plasma medium [5]. 

It is a well known fact that cold atomic gases confined 
in MOTs sustain various types of instabilities. One ex- 
ample is found in the work realized by Kim et al. [9], 
where a parametric instability is excited by an intensity 
modulated laser beam. This instability is based on a 
nonlinear response of the individual atoms to one of the 
control parameters of the MOT. On the other side, spon- 
taneous instabilities arising for large MOTs have been 
described in p!0Hl2] , where the feedback of retroreflected 
laser beams of the MOT can induce stochastic or deter- 
ministic chaos for a large optical thickness of the MOT. 
These instabilities have been investigated by looking at 
the motion of the center of mass of the cloud. A further 
type of spontaneous instabilities, due to the large spatial 
extend of the MOT has been described in [16]. Here, the 
instability finds its origin at the edge of the MOT where 
the Zecman shift produces a local change of the detun- 
ing seen by the atoms with a negative friction for atoms 
beyond a certain horizon. 

In the present work, wc theoretically explore the possi- 
ble existence of a new instability process, due to the cou- 
pling between two collective oscillations in a MOT: the 
"plasma" and the center-of-mass (c.o.m.) modes. We use 



the formalism developed in our recent work [S] , where we 
have shown that plasma-hybrid waves can propagate in 
the atomic cloud and have established their dispersion re- 
lation. These hybrid waves exhibit features that are com- 
mon to both electron-plasma and acoustic waves [T3HT5] . 
Here we start with a set of fluid equations describing the 
dynamics of the atomic medium, where we have retained 
the main features of the cooling and trapping forces. We 
then show that these plasma-like oscillations can couple 
with the center-of-mass of the cloud, giving place to the a 
possible parametric instability mechanism. By identify- 
ing the relevant parameters and establishing the stability 
criteria, we are able to show that this mechanism could 
lead to an instability threshold in the same parameter 
range as the one observed by Labeyrie et al. [16]. New 
instability regimes are also predicted. 

Following our previous work ^ , the collective behavior 
in the trap is governed by a set of fluid equations 
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where n and v represent the gas density and velocity, 
respectively, and M is the atomic mass. The total force 
is given by F = Fc -I- Ft , where Fc is the collective force 
and Q — {cr r — a L)cf lIq / c represents the square of the 
effective electric charge [71 [8] , with c being the speed of 
light, /g the total intensity of the six laser beams and uji 
and (Tl representing the emission and absorption cross 
sections [3J. The term Ft — — VJ7 stands fot the trapping 
force. We consider a low intensity Doppler regime MOT, 
where the trapping potential is given by 
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Here, fiB represents the Bohr magneton, k = k{S,Io/Is) 
is the spring constant, a = a{S, lo/h) is the friction coef- 
ficient, 6 is the laser detuning, 1^ is the atomic saturation 
intensity and VB = |Vi3| represents the magnetic field 
gradient [5]. In such a picture, the gas confined in a MOT 
can be regarded as an one-component plasma, where the 
effective plasma frequency is given by 



(5) 



It is assumed that the equilibrium density of the gas uq = 
no{r) depends only the radial coordinate r and angular 
distortions of the cloud will be ignored. For simplicity, 
we will use a water-bag model, where we assume that the 
density is approximately constant inside the cloud. 

We linearize Eqs. ([l])-([3]) by separating each relevant 
physical quantity into its equilibrium and perturbation 
components, such that 



71 = no +711+ 7^2, V = Vi+V2, F = Fi+F2. (6) 

The subscripts 1 and 2 label two different perturbations, 
corresponding to two distinct collective modes in the 
MOT. This procedure allows one to include any pair of 
modes and could be extended to more than two modes, 
thus broadening the range of applicability of the model. 
In this work, we consider that 1 and 2 label the per- 
turbations due to center-of-mass (c.o.m.) and plasma 
oscillations. Accordingly, Fi = Fj is the restoring force 
in the c.o.m. and F2 = Fc represents the collective force 
responsible for the plasma oscillations. For moderate os- 
cillation amplitudes, the centcr-of-mass velocity Vi is de- 
scribed by a harmonic oscillation 
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+ WcAf"-2 COsilVcAlt + 4>)] = 0. 



The third term of this expression describes the nonlin- 
ear coupling between the two collective modes mentioned 
above. This is one of the main results of the paper. No- 
tice that the plasma mode ^ could be replaced by any 
of those in the hierarchy w„ ; associated to the Tonks- 
Dattner (TD) resonances, the standing waves inside the 
cloud which hierarchy provides a large number of modes 
with frequencies very close to each other [8]. The reason 
for this choice is because wp is the lowest mode in the 
TD series. 

Making use of the ansatz ?i2(r,t) — B2{r)A2{t), Eq. 
^ yields 



[v + 2e cos(2t)] ^2 + 6 sin(2r' 
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where ^ = 2ui • V In ,62 is a free nonlinear coupling pa- 
rameter and describes the amplitude of the center-of- 
mass oscillation. Here, 2t = wcm^ + ^ — '^^p/^cm 
and e — ^/wcm are dimensionless variables. Equation 



( 10 1 describes the dynamics of a parametrically excited 
system and belongs to the family of Hill equations. It is 
formally similar to the Mathieu equation, which is well- 
known in the literature for containing unstable solutions. 
According to the standard Floquet theory 23 , the solu- 
tions to Eq. ( 10 ) can be expanded into Fourier series 
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vi(i) = ui sin{oJcMt + ' 
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where ojcm (associated to Fi - for the sake of generality, 
the details on the calculation of locm will be discussed 
later on) represents the c.o.m. mode frequency, Ui its 
amplitude and is an arbitrary phase. We will consider 
that the c.o.m. oscillates as a rigid body, which means 
that the perturbation 711 can be neglected (iii <C 712 <C 
7io). Combining these approximations with the above 
fluid equations, one easily obtains 
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+ n2UJcM cos{ujcMt + <P)] = 0, 
Noticing that Tio^tV • V2 = ujpn2, it follows 
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where 7 is the characteristic exponent |23j . Plugging it 
into Eq. (10), one obtains a system of equations for the 
coefhcients a„. The non-trivial solution is obtained by 
computing the roots of the so-called infinite Hill deter- 
minant, i.e.. 



det(Aj-fc(7i)) = det{[i^ + {j + 2jn)^]Sj.k 



(12) 
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where Sj,k is the Kroenecker delta. In order to carry out 
calculations, only the truncated system 71 — ~N,...,N 
is considered. This approximation results on the Ince- 
Strutt diagram jM]. In Fig. ([T]) we plot the stability 
chart of Eq. ( 10 ) in the — e plan, truncated at iV = 5, 



where the vr and 27r-periodic marginal curves correspond 
to 7 = and 7 = 7, respectively. The marginal curves 
i/(e) separate the different regions of stability (check Ref. 
P5] for further details on the Floquet theory). 



3 



For comparison with experiments, it is necessary to 
express the frequencies uj p and uc m in terms of real-Ufe 
parameters. The balance between the trapping and col- 
lective forces is given by the hydrodynamical equilibrium 
condition D/Dtv = {d/dt + v • V)v — 0, which simply 
corresponds to set 



0, 



(13) 



where Fj = — V[/ = —ky. Taking the divergence of the 
later, we obtain the following relation 
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which simply corresponds to the expression for the com- 
pression limit of the MOT [5]. This condition estab- 
lishes a relation between the plasma ([s]) and the c.o.m 
oscillation frequencies: ujcm = i^t = i^p/V^, where 
(jJt = \J njM is the trapping frequency. The identity 
'^CA/ — '^t is equivalent of setting Fi = Ft and it is 
simply a consequence of the Kohn's theorem (17j , which 
basically states that interactions do not affect the c.o.m. 
motion. The equality IiJcm — wp/^/i is often recognized 
in unscreened Coulomb systems [THl [12] , as it describes 
the so-called Mie resonance. The analogy between MOTs 
and Coulomb systems follows from the fact that the effect 
of the trap is replaced by the background ionic density in 
the later. For typical experimental conditions of a cold 
^^Rb gas, operating at 5 = — l.SF, where F 6 MHz 
is the natural lifetime of the atomic transition, with a 
magnetic field gradient VB = 5 G/cm and Iq/Is ~ 0.3, 
the corresponding plasma frequency is ujp/2tt « 120 Hz. 
As a consequence of the condition (14), we set the value 
V ^ 12 in Eq. (10) and, correspondingly, the instability 



occurs provided that 



ei < e < £2 and e > £3, 



(15) 



where ei = 3.17, £2 = 3.91 and £3 = 4.91 (see Fig. ([T])). 
Using the definition of £, and restricting the discussion to 
the low saturation limit, the marginal curves f = ^(£j)i 
with i = (1,2,3), representing the critical values of the 
nonlinear coupling parameter, are simply given by 
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where fc^ is the laser wave vector. In Fig. ([2| we plot the 
marginal curves ^1, ^2 and ^3 against the relevant MOT 
parameters. Below the marginal curve ^1, the oscilla- 
tions are stable. However, for ^1 < ^ < ^2, the system 
undergoes large amplitude density oscillations and the 
becomes unstable. The oscillations become stable be- 
tween the marginal curves ^2 and ^3. Finally for ^ > ^3 
we can observe that the oscillations are again unstable. 



This way of coupling different modes of the system, 
added to an external drive (or sufficient power at the 
correct frequency in the noise spectrum) and therefore 
controlling the value could lead to an instability based 
on mechanisms similar to the ones used in [5]. Applying 
this idea to a typical range of {5, VB) in Eq. (16), it 



is possible to build a stability chart similar to that pre- 
sented in Fig. We note that a stability diagram 
similar to the one observed by Labeyrie et al. |16j can be 
obtained for specific values of driving parameters, corre- 
sponding here to ^ = 280 Hz. In the present paper, how- 
ever, the basic ingredient for the instability is different 
from that used in [16], what therefore suggests it is pos- 
sible to distinguish these two mechanisms: in the later, 
the instability is due to a change of sign in the friction 
coefhcient at the edge of the cloud [52] , where the driv- 
ing motion takes place. Instead, we here suggest that an 
instability in the same parameter range can be obtained 
by coupling two (or more) collective modes, where the 
energy is transfered from one to the other, as an example 
of a density-wave mixing phenomenon. Our model also 
predicts different stability zones, which does not happen 
in the model of Ref. [16 . We believe that the main dif- 
ference is the fact that our model does not depend on 
the size of the cloud, which is consistent with the fluid 
approximation. As a consequence, we do not require the 
edge to contribute to the dynamics of the system. We 
also argue that the coupling between the c.o.m and the 
different TD modes may occur in real experiments (as 
they are very close to ujp) and it should be responsible 
for the broadening of the intermediate instability zone 
Ui in Fig. [3j which better agrees with the experimental 
observations of Ref. [T^. This is because for ly > 12 the 
intermediate instability zone £1 < £ < £2 is larger, as one 
can directly observe in Fig. ([T]). 

A short discussion of the validity range of the present 
model is in order. We expect the present results to 
be valid only valid in the low-intensity Doppler regime, 
where the dynamics of the single atom MOT remains es- 
sentially linear. The existence of nonlinear corrections is 
obvious from the above fluid equations, which are not at 
all related with the nonlinearity in the basic forces, but 
are only rooted in the collective behavior of the MOT. 

In this work, we have highlighted several aspects of 
the fluid description of the collective behavior of a cold 
atomic gas in a magneto-optical trap. We have estab- 
lished the governing equation for the nonlinearly coupled 
dynamics of the center-of-mass and the plasma oscilla- 
tions of the cloud, which we have shown to be the root 
of an instability mechanism. It was also shown that the 
condition of dynamical equilibrium provides a relation 
between the frequencies of the center-of-mass and the 
plasma oscillations, ujcm = uip/V^. Our model yields 
an instability threshold which can be close to the one 
observed in [16j . even though mechanisms not included 
in this work have been put forward to explain these 
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FIG. 1: Stability chart of Eq.(llO|, obtained for a n = 5 Hill 
determinant. The full lines represent 7r-period solutions and 
the dashed lines represent the 27r-period ones. The shadowed 
(light) regions correspond to stable (unstable) solutions. The 
red dots represent the marginal values ei , 62 and es discussed 
in the text. 
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FIG. 2: (Color online) Marginal curves = ^(e^) (i = 1, 2, 3) 
as function of: a) the detuning 5, at V-B = 3 G/cm and 
b) the magnetic field gradient Vi?, sX 6 = — O.TSF. In both 
cases, lo/Is = 0.3 and (solid line), ^2 (dashed line) and ^3 
(dotted line) are the marginal curves discussed in the text. 
The shadowed areas are stable. 



self-sustained instabilities. In addition, we predict the 
existence of new stability regions, compatible with the 
same coupling parameter, thus giving a new insight to 
the problem and proposing new areas of experimental 
research. We concluded that our model can be used to 
describe, at least qualitatively, a family of instabilities of 
the collective oscillations in a MOT. 
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